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Joint  Agency  Turbulence 
Experiment  Interim  Report 


1.  INTRODUCTION 

The  Air  Force  Geophysics  Laboratory  (AFGL)  is  developing  an  airborne 
remote  turbulence  sensor  to  enable  pilots  to  avoid  regions  of  turbulence  that  are 
hazardous  to  aircraft.  Methodologies  developed  at  AFGL  were  initially  tested 
with  ground-based  radar  and  aircraft  gust  data  acquired  from  the  1976  Rough  Rider 
program.  During  this  experiment,  the  National  Severe  Storms  Laboratory  Norman 
Doppler  radar  performed  periodic  sector  scans  of  a  storm  while  an  Air  Force  F4 
Phantom  aircraft,  equipped  with  lightning  and  vertical  gust  measurement  equip¬ 
ment,  made  storm  penetrations.  During  these  penetrations  aircraft  measurements 
were  made  from  which  the  vertical  gust  velocities  were  derived. 

These  radar  and  aircraft  data  were  previously  analyzed  and  reported 

12 

(Bohne  '  ).  The  results  of  these  observations  showed  definite  promise  in  the  use 
of  radar  to  detect  and  quantify  regions  of  turbulence  in  storms. 


(Received  for  publication  1  July  1983) 

1.  Bohne,  A.  R.  (1981a)  Radar  Detection  of  Turbulence  in  Thunderstorms. 

AFGL-TR-81-0102,  ADA108679. 

Bohne,  A.  R.  (1981b)  estimation  of  turbulence  severity  in  precipitation 
environments  by  radar.  Preprints,  20th  Radar  Meteorology  Conference 
Amer.  Meteor.  Soc.  ,  Boston,  Mass.,  pp.  446-453. 


To  increase  the  data  base  and  resolve  some  of  the  discrepancies  observed  in 
this  early  data,  a  more  coordinated  set  of  data  was  required.  To  this  aim  AFGL 
established  the  Joint  Agency  Turbulence  Experiment  and  participated  in  the  Storm 
Hazards  Program  at  Wallops  Island,  Virginia,  with  NASA  Langley  Research  Center 
and  the  NASA  Wallops  Flight  Center  during  the  early  summer  months  of  1981  and 
1982.  This  effort  utilized  ground-based  and  tracking  radars,  a  lightning  and  gust 
measurement  equipped  F106B  aircraft,  and  other  ground-based  sensors  to  study 
the  effects  of  lightning  and  turbulence  on  aircraft. 

The  particular  aim  of  AFGL  was  to  obtain  ground-based  radar  and  in-situ 
aircraft  gust  data  that  were  spatially  and  temporally  coincident.  The  AFGL 
incorporated  Doppler  processing  and  recording  equipment  into  the  Wallops  Island 
SPANDAR  radar  system.  This  radar,  with  a  0.41°  full  half-power  beam  width  and 
wavelength  of  10.  56  cm,  was  operated  at  320,  640,  or  1280  prf,  which  thus 
afforded  unambiguous  ranges  and  velocities  of  468,  234,  and  117  km  and  33.  79. 
16.90,  and  8.  45  m/s,  respectively. 

The  time  and  space  coincident  aircraft  and  radar  data  acquired  during  these 
tests  are  being  used  to  further  test  the  AFGL  turbulence  detection  methodology, 
and  to  provide  useful  information  on  the  state  of  the  turbulent  atmosphere.  The 
data  presented  here  are  observations  of  the  data  structure  and  environments  Irom 
the  1981  season.  The  integrated  radar-aircraft  analysis  results  will  appear  in  a 
later  report. 

2.  OPERATIONAL  MODE 

Two  different  radar  data  sets,  tracking  and  general  coverage,  were  recorded 
by  AFGL.  During  aircraft  storm  penetrations,  the  aircraft  was  tracked  by  the 
SPANDAR  radar.  This  was  accomplished  by  slaving  SPANDAR,  parallax  cor¬ 
rected,  to  a  nearby  tracking  radar.  Using  the  aircraft  position  information  sup¬ 
plied  by  the  tracking  radar,  six  contiguous  tracking  gates  were  locked  onto  the 
aircraft  position.  The  aircraft  was  kept  in  a  selected  tracking  gate  to  ensure  that 
the  aircraft  track  was  being  maintained.  The  tracking  gate  data  recorded  were 
the  pulse-to-pulse  in-phase  and  quadrature  Doppler  returns  and  ancillary  informa¬ 
tion.  Real-time  monitoring  was  accomplished  by  sampling  tracking  gate  data  \  i t h 
a  Ubiquitous  Spectrum  Analyzer  and  observing  the  aircraft  spike  in  the  one- 
dimensional  Doppler  spectrum.  The  data  were  recorded  in  analog  mode  on  an 
Ampex  FR1300  recorder. 

Between  storm  penetrations,  AFGL  collected  Doppler  radar  lata  in  sector 
scan  mode.  Typically  768  or  1024  1-psec  gates,  contiguous  or  with  a  l-p.-e-c 
spacing,  were  scanned  over  a  30°  to  60°  sector  with  three  to  five  elevation  steps. 


The  reflectivity  factor  and  Doppler  spectrum  mean  and  variance  were  developed 
in  real  time  by  the  hardwired  pulse-pair  processor  and  recorded  on  the  AFGL 
archiving  system.  The  sector  scan  data  were  used  to  provide  real-time  aircraft 
hazard  avoidance  information  to  aid  in  selecting  the  penetration  routes  for  the 
F106B,  and  to  provide  general  storm  surveillance. 

During  storm  penetrations  the  F10GB  equipment,  in  addition  to  making  exten¬ 
sive  lightning  strike,  atmospheric  chemical,  and  multi-wavelength  observations, 
made  gust  parameter  measurements  of  the  three  components  of  wind  gust  velocity. 
These  measurements,  made  at  the  rate  of  80/s,  were  later  sampled  and  digitized 
at  the  rate  of  20/s.  The  measurements  were  then  used  to  derive  the  longitudir 
transverse,  and  vertical  gust  velocity  estimates.  These  data  will  provide  the 
'ground  truth'  measurements  against  which  the  radar  data  will  be  compared,  ai 
will  provide  further  information  on  the  state  of  the  atmosphere,  which  is  diffic 
to  assess  from  the  radar  data  alone. 

3.  DATA  REDUCTION 

The  tracking  gate  data  consisted  of  pulse-to-pulse  measurements  of  the  in- 
phase  and  quadrature  sampled  return  and  associated  housekeeping.  These  track¬ 
ing  data  were  recorded  in  analog  fashion  on  fourteen  track  tape.  Ail  data  wore 
initially  reviewed  with  use  of  a  Ubiquitous  Spectrum  Analyzer  to  identify  those 
gates  in  which  aircraft  contamination  was  present.  One  of  the  gates  from  each 
penetration  was  then  digitized  for  use  in  later  analyses.  Typically,  the  uncon¬ 
taminated  gate  between  the  radar  and  the  aircraft  closest  to  the  plane  was  used. 

It  was  observed  that,  with  the  aircraft  located  near  the  middle  of  the  six  contig¬ 
uous  tracking  gates,  the  first  and  last  pair  of  gates  (1,2  and  5,  (5)  were  not  con¬ 
taminated.  The  selected  data  were  then  digitized  and  recorded  as  25fi  (128  real 
and  imaginary)  sample  records.  For  comparison  purposes,  three  different  esti¬ 
mators  were  employed  to  determine  Doppler  spectrum  mean  and  variance.  These 
included  pulse  pair,  objective  thresholding,  and  constant  threshold  from  the 
spectrum  peak.  The  algorithms  employ  standard  autocorrelation  and  FFT  tech¬ 
niques  and  include  automatic  image  spectrum  removal  and  velocity  unfolding. 

Because  the  aircraft  was  frequently  in  precipitation  sparse  regions,  small 
segments  of  missing  data  were  occasionally  observed  in  the  time  histories  of  the 
track  data.  Analyses  of  the  radar  radial  velocity  data  include  autocorrelation, 
structure  function,  and  power  spectrum  analysis  methods.  To  implement  these 
most  easily  it  is  necessary  to  have  continuous  time  histories  of  data.  To  remove 
the  gaps  (one  to  six  missing  consecutive  data  values)  the  mean  value  and  variance 
of  the  parameter  in  question  was  estimated  for-  a  small  neigiibv  food  of  values  on 


4 

i 

* 


i 


each  side  of  the  data  gap.  The  mean  values  established  a  linear  trend  acr  oss  the 
data  gap.  The  trend  values  were  then  modified  by  adding  a  random  value  at  each 
missing  data  location.  The  family  of  random  values  used  had  a  variance  that  was 
similar  to  that  observed  in  the  neighboring  areas.  The  fitted  r  egions  were  found 
visually  undetectable  against  the  background  of  true  neighboring  data.  Although 
this  technique  does  not  replicate  the  true  missing  data,  the  fitted  segments  are  so 
small  in  number  and  extent  in  contrast  to  the  original  total  time  histories  that 
little  distortion  of  energy  content  at  the  frequencies  of  concern  are  expected. 

Further  editing  was  occasionally  required  to  remove  spurious  spectral 
estimates  believed  to  result  from  occasional  wavering  of  the  tracking  gate  location 
about  the  aircraft  position.  The  AFGL  system  used  the  tracking  information 
supplied  to  the  SPANDAR  radar  to  update  the  tracking  gates  at  a  rate  of  10  s. 

This  time  interval  represents  an  aircraft  movement  of  roughly  10  m  between  up¬ 
dates.  This  is  small,  however,  in  comparison  to  the  AFGL  tracking  range  resolu¬ 
tion  of  60  m. 

Of  added  concern  are  the  effects  of  the  occasional  large  ranges  of  the  storms 
penetrated.  Here,  small  pointing  errors  could  easily  move  the  tracking  gates  by 
a  few  hundreds  of  meters.  The  combination  of  the  range  stepping  update  of  the 
AFGL  system  with  a  slight  pointing  error  could  result  in  the  time  series  data 
being  used  to  derive  a  single  Doppler  spectrum  to  be  obtained  from  nearly  sepa¬ 
rate  spatial  locations,  possibly  exhibiting  different  wind  characteristics.  This  is 
believed  to  be  primarily  responsible  for  the  spurious  results  observed.  The 
above-mentioned  procedures  were  performed  on  both  the  Doppler  spectrum 
velocity  and  variance  data. 

4.  DATA  PRESENTATION 
4.1  Sounding  Data 

The  environmental  soundings  to  be  shown  represent  estimr-'es  of  die  structure 


at  the  time  and  location  of  observations.  Generali-.,  'he  ohservn'i  >r.  n.mio  Is 
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4.2  Cartesian  Representation  of  Sector  Scan  Data 

The  data  acquired  during  the  sector  scan  portion  of  storm  observations  are 
interpolated  to  a  Cartesian  grid  system.  The  data  are  the  equivalent  reflectivity 
factor  and  the  Doppler  mean  velocity  at  storm  penetration  altitude.  This  level 
was  determined  from  the  radar  computed  tracking  altitude,  rather  than  from  the 
aircraft  determined  altitude,  to  ensure  proper  accounting  for  Earth  curvature  and 
index  of  refraction  effects.  The  reflectivity  factor  and  radial  velocity  structures 
were  observed  over  a  few  sector  scans  where  possible  to  determine  an  advection 
velocity  for  each  storm  during  the  observation  period.  This  advection  velocity 
was  then  applied  to  all  scan  data  to  ensure  that  the  data  were  interpolated  to  the 
proper  grid-point  locations.  In  all  cases  this  advection  velocity  is  in  good  agree¬ 
ment  with  the  estimated  sounding  velocity  for  the  observational  period. 

Often  a  number  of  aircraft  storm  penetrations  were  performed  before  a  com¬ 
plete  sector  scan  update  could  be  obtained.  In  this  instance,  a  single  represen¬ 
tative  Cartesian  data  set,  with  the  aircraft  penetration  tracks  advected  to  the 
proper  storm  locations  as  determined  by  the  storm  advection  velocity,  is  used. 
Now,  however,  comparisons  between  aircraft  and  radar  scan  data  do  require  a 
relatively  stationary  storm  structure.  Those  cases  where  significant  evolution 
was  observed  between  the  time  of  aircraft  storm  penetrations  and  the  sector  scan 
period  are  noted.  The  grid-point  separations  in  these  plots  range  from  0.  5  to 
1.  5  km  in  the  horizontal  and  0.  5  to  1.0  km  in  the  vertical.  The  occasional  coarse 
resolution  is  a  result  of  the  limited  number  of  sector  scans  performed  between 
penetrations  and  the  large  range  of  the  storms  from  the  SPANDAR  radar. 

Finally,  since  a  15-dBZ  threshold  is  used  in  the  plotting  of  the  reflectivity  data, 
the  plots  do  not  represent  the  entire  spatial  extent  of  the  storms. 

4.3  Tracking  Data 

The  tracking  gate  Doppler  data  displayed  are  the  Doppler  spectrum  mean 
velocity  and  variance.  The  data  presented  are  a  smoothed  version  of  the  original 
data  and  were  obtained  by  applying  a  simple  constant  nine-point  filter  on  the  nine 
nearest  neighbor  values  centered  about  the  data  point  in  question.  The  smoothed 
data  show  more  clearly  the  trend  in  the  data  and  allow  easier  comparison  of  radar 
spectrum  data  with  storm  structure  encountered  along  the  aircraft  track.  This 
filter  will  not  be  employed  in  future  power  spectrum  analysis. 

4.4  Turbulence  Severity  Estimates 

The  parameter  of  interest  here  is  the  eddy  dissipation  rate  (e),  which  repre¬ 
sents  the  rate  of  transfer  of  turbulence  kinetic  energy  from  larger  to  smaller 


scales  of  motion,  or  eddies.  The  turbulence  severity  index,  simply  the  cube  root 
of  e,  has  been  offered  (MacCready  )  as  a  means  of  quantifying  the  degree  of 
hazard  to  aircraft  penetrating  a  given  region  of  turbulence.  Radar  estimates  of 
eddy  dissipation  rate  (e)  along  the  aircraft  penetration  tracks  are  determined  us¬ 
ing  two  different  methods.  First,  the  tracking  gate  Doppler  velocity  data  are 
used  in  a  structure  function  analysis,  which  essentially  measures  the  degree  of 
correlation  of  radar  radial  velocity  data  with  distance.  Second,  the  tracking  gate 
Doppler  spectrum  variance  (DSV)  data  are  related  to  the  eddy  dissipation  rate  via 
a  relation  founded  on  the  concept  that  the  turbulent  air  and  precipitation  motions, 
particularly  at  scales  primarily  less  than  the  maximum  pulse  volume  dimension, 
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are  mapped  into  DSV.  The  relation  used  is  that  derived  by  Frisch  and  Clifford, 
which  assumes  a  turbulence  outer  scale  (effective  maximum  eddy  size  in  the 
inertial  subrange)  of  infinity.  The  more  accurate  method  which  employs  use  of 
actual  estimates  of  the  turbulence  outer  scale  to  relate  spectrum  variance  and  e 
were  not  made  at  this  time.  The  following  material  briefly  outlines  the  manner  in 
which  the  two  methods  were  employed. 

The  form  of  the  one-dimensional  structure  function  used  is 

D(r)  =  (D^r)  -  Dt(r))(xf/r)2  +  Dt(r)  '  (1) 

=  <(V(R)  -  V(R  +  r))2)  .  (2) 

where  D^  (r)  and  Dt(r)  are  the  one-dimensional  longitudinal  and  transverse  struc¬ 
ture  functions,  respectively,  and  are  given  by 

D(|r)  =  C(fr)(2/3) 4  (3) 

Dt(r)  =  4/3  C(e  r)(2/3)  ,  (4) 

for  the  inertial  subrange  region  of  a  homogeneous  isotropic  turbulence  field.  This 
relation  is  required  since  the  aircraft  penetration  tracks  were  neither  parallel 
nor  orthogonal  to  the  radar  viewing  direction.  In  these  relations  C  is  a  universal 
constant  (1.77),  r  is  the  distance  interval  between  successive  radar  velocity 
measurements  used  to  form  the  estimate,  and  is  the  component  of  this  distance 
along  the  radar  viewing  direction. 

3.  MacCreadv,  P.  (1964)  Standa rdization  of  gustiness  values  fr  om  aircraft, 

J.  Appl.  Meteorol,  _13:808-811. 

4.  Frisch,  A.  5.,  and  Clifford,  S.  P.  (1974)  A  study  of  convection  capped  by  a 
stable  Inver  using  Doppler  radar  and  acoustic  echo  sounders,  J.  Atmos,  Sri 
31:1922-1928. 
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when  using  ra-inr  ■  ■ ; •  1  i n  1  velocity  values,  the  random  fluctuation  in  Donnie:-  >■  i  >. 
arising  front  turbulence  scales  less  than  the  maximum  pulse  voluit.  -  ’.intension 
may  be  strongly  filtered  out.  Thus,  the  radar  structure  function  estimates, 
least  for  the  case  of  an  ensemble  of  measurements,  will  be  diminished  in  magni¬ 
tude  and  contain  proportionately  greater  amounts  of  energy  from  the  larger  scales 
of  motion. 

The  mode  of  movement  of  the  tracking  gates  alternated  from  one  of  nearly 
uniform  scanning  to  a  1-psec  range  jump  as  the  plane  moved  in  or  out  in  range. 
Thus,  the  tracking  gate  data  comprising  the  continuous  time  histories  are  not 
uniformly  spaced.  To  obtain  a  more  spatially  uniform  data  set  for  use  in  the 
structure  function  analyses,  only  those  data  that  were  obtained  immediately  after 
a  range  jump  are  employed.  This  results  in  basic  data  sets  where  successive 
velocity  values  are  150  to  240  m  apart.  The  structure  function  estimates  at  each 
data  point  location  are  derived  from  a  local  data  set  centered  around  the  data 
point  location  of  interest.  These  local  data  sets  varied  in  length  from  3  to  5  km. 

The  estimates  of  e  determined  by  this  method  will  exhibit  some  error.  First, 
the  structure  function  estimates  are  derived  from  a  limited  number  of  observa¬ 
tions  that  may  only  loosely  approximate  an  ensemble  of  observations.  Second,  the 
magnitude  of  e  is  diminished  due  to  pulse  volume  filtering  action.  Third,  the 
estimates  are  biased  through  incorporation  of  nonturbulent  storm  stiucture  wind 
shear.  Fourth,  each  velocity  value  naturally  includes  a  random  error  component 

related  to  Doppler  spectrum  breadth  and  sample  dwell  time.  This  may  add  from 
2 

0.  1  to  0.4  (m/s)  to  the  structure  function  magnitude. 

In  an  attempt  to  account  for  these  effects,  the  local  data  sets  are  first  linearly 

detrended  and  eddy  dissipation  rate  estimates  are  determined  at  the  smallest 

2  /3 

length  scales  where  the  r  behavior  is  reasonably  well  observed,  typically  at 
length  scales  of  300  through  500  m. 

The  second  method  employed  to  estimate  e  relies  upon  the  relationship 
between  DSV  and  the  turbulent  air  motions.  This  method  assumes  some  knowledge 
of  the  effective  size  of  the  largest  turbulence  eddies  within  the  inertial  subrange. 
The  term  effective  is  used  since  we  are  replacing  the  actual  turbulence  energy 
spectrum  that  contains  energy -containing  eddies,  in  addition  to  the  inertial  sub¬ 
range,  with  a  fully  Kolmogorov  energy  spectrum  form.  If  the  actual  Von  Karman 
form  of  energy  spectrum  is  known,  or  fine  scale  in-situ  measurements  are  avail¬ 
able,  then  the  effective  outer  scale  length  may  be  determined.  These  data  are  not 
presently  available,  however.  Thus  the  solution  of  Frisch  and  Clifford,  ^  whereby 
the  inertial  subrange  is  assumed  to  extend  to  eddies  of  infinite  size,  is  employed. 
The  relationship  used  is 
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Var  =  1.354  C  (e  a.V  1  '  (1  -  z/15  -  (z“)/105  -  ...) 


(5) 


where  C  is  a  universal  constant  of  value  near  1.35,  and  z  and  a  are  related  to 
effective  pulse  volume  width  (R0)  and  length  (L)  by 


z  =  1  -  (L/R0)2 

(fi) 

a  =  R0/(8  In  4)(1/2)  . 

(7) 

At  the  large  ranges  investigated,  the  series  term 
0.93,  yielding  the  working  relation 

remains  essentially  constant  at 

Var  =  .0225  C  <e  R>(2/3)  . 

(8) 

The  variance  data  used  in  these  analyses  are  the  companion  data  to  the 
Doppler  velocity  data  used  in  the  structure  function  analyses.  The  DSV  values 
are  derived  from  block  averaging  all  variance  data  contained  in  five  successive 
range  jump  segments.  Effectively,  this  is  simply  averaging  DSV  over  a  spatial 
range  comparable  to  that  used  in  the  structure  function  analyses.  Eddy  dissipation 
rate  is  then  determined  using  these  averaged  variance  values. 

It  is  important  to  note  that  the  structure  function  and  spectrum  variance  tech¬ 
niques  rely  upon  different  effects  of  the  turbulent  air  (precipitation)  motions  on  the 
radar  Doppler  spectrum  parameters.  First,  the  structure  function  relies  upon 
the  turbulent  air  motions  being  mapped  into  fluctuation  of  the  Doppler  velocity,  and 
then  relates  this  fluctuation  with  intensity  of  the  turbulence  field.  In  the  form  used 
here  [Eq.  (1)],  no  accounting  for  the  pulse  volume  filtering  effects,  that  is.  trans¬ 
ference  of  some  of  this  fluctuation  energy  into  Doppler  spectrum  breadth  is  per¬ 
formed,  The  variance  method,  on  the  other  hand,  relies  upon  the  spectrum 
broadening  effects,  and  does  account  for  loss  of  measurable  turbulence  energy  to 
the  fluctuating  mean  velocity.  However,  the  relation  used  here  does  not  assume 
a  finite  outer  scale  length.  Past  experience  has  shown  that  a  length  A  of  0.  5  to 
1. 5  km  is  usually  appropriate,  indicating  use  of  Eq.  (5)  will  underestimate  f .  If 
the  outer  scale  length  is  considerably  lar  ger  than  twice  the  maximum  pulse  volume 
dimension,  then  Eq.  (5)  is  quite  good. 

Neglecting  for  the  moment  the  various  beam  filtering,  storm  structure  wind 
shear,  and  other  bias  and  random  errors  discussed  earlier,  comparison  of  the 
time  history  patterns  and  the  relative  magnitudes  of  the  two  separate  estimates 
of  e  may  provide  some  useful  qualitative  information  on  the  turbulence  sever  ity, 
turbulence  outer  scale  length,  and  presence  of  strong  storm  wind  shear'. 


Hi 


f .  1 A.V.V^iO f «•. 


For  example,  the  maximum  pulse  volume  dimension  is  typically  0.  6  to  1,  4  km. 
With  an  ensemble  of  measurements,  about  90  percent  of  the  turbulence  energy  is 
mapped  into  DSV  for  turbulence  fields  having  effective  outer  scale  lengths  up 
to  2.  8  km.  This  is  certainly  within  the  range  of  expected  values.  The  structure 
function  estimates  employ  velocity  values  that  are  spatial  averages  over  the  pulse 
volume  (roughly  150  by  1  km),  and  which  are  separated  in  distance  by  150  to 
240  m.  If  the  turbulence  is  homogeneous  and  isotropic,  and  the  structure  function 
estimates  (reasonably  well)  approximate  an  ensemble  of  measurements,  then  we 
would  expect  most  turbulence  energy  to  be  mapped  into  DSV,  with  little  mapped 
into  the  fluctuation  of  the  mean  velocity.  Thus,  if  the  meteorological  field  is 
basically  turbulent  and  the  outer  scale  length  less  than  twice  the  maximum  pulse 
volume  dimension,  then  we  expect  the  time  histories  of  the  two  estimates  to  be 
similar  in  form,  but  with  the  structure  function  (STF)  derived  estimates  to  be 
smaller  in  magnitude  than  the  DSV  values.  If  the  magnitudes  are  small  and  equal, 
and  the  patterns  are  similar,  then  the  effective  outer  scale  may  be  small  enough 
to  allow  the  DSV  underestimate  to  approach  the  STF  value. 

If  the  DSV  estimate  varies  in  magnitude  while  the  STF  estimate  is  constant, 
then  a  change  in  the  turbulence  outer  scale  is  indicated.  If  the  DSV  estimate  is 
relatively  constant,  while  the  STF  estimate  varies  significantly,  then  presence  of 
wind  shear  biasing  the  STF  estimate  is  suggested.  This  results  from  the  fact  that 
a  simple  transverse  wind  shear  contribution  to  e  is  approximately  eight  times 
greater  for  the  STF  than  the  DSV  estimate  at  these  ranges. 

Thus,  with  this  set  of  qualitative  guidelines,  and  remaining  mindful  of  the 
caveats  mentioned  earlier,  some  useful  information  concerning  the  turbulence 
structure  may  be  derived.  Finally,  one  must  note  that  turbulence  severity 
(c^/3))  j^s  displayed,  and  differences  in  the  magnitudes  of  e  from  the  STF  and 
DSV  estimates  are  much  greater  than  shown  on  these  plots. 

5.  DISCUSSION  OF  DATA 
5.1  1  July  1981 

The  first  data  set  was  obtained  on  1  July  1981  (day  182).  Figure  1  shows  the 
estimated  environmental  sounding  data.  The  wind  direction  varies  slowly  from 
about  170°  near  the  surface,  to  120°  at  (5  km.  Over  this  height  the  speed  also 
increases  slowly  from  about  14  to  20  m  s.  The  temperature  data  is  suspect  since 
the  sounding  shows  the  environment  becoming  very  dry  at  3 -km  height,  whereas 
the  penetration  level,  as  determined  by  the  F10BB,  was  4.  57  km. 

Figure  2  shows  the  reflectivity  factor  for  the  storm  at  the  penetration  altitude 
of  the  F106B.  The  time  period  over  which  the  sector  scan  data  were  collected 
extends  from  18:18:00  to  18: 17:28  GMT.  The  plot  is  time-adverted  to  repr  esent 
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Figure  1.  Wallops  Island  Sounding  Data  for  1  July  1981  at  1200  GMT 
Wind  data;  speed  (dash)  and  direction  (solid):  temperature;  tempera¬ 
ture  (solid),  dew  point  (dash),  and  pseudoadiabat  (dash-dot) 
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Figure  2.  Contours  of  Storm  Reflectivity  Factor  (dBZ)  for  1  July  1981  on 
Constant  Height  Surface  at  Aircraft  Storm  Penetration  Altitude  of 
4.  57  km.  Contours  are  in  5-dBZ  increments  with  a  recycling  dot,  dash, 
and  solid  pattern.  Minimum  contour  value  is  10  dBZ.  Time-adjusted 
aircraft  track  is  shown.  Distances  are  relative  to  SPAMDAR  radar. 

Time  is  18:17:00  GMT 


the  expected  appearance  at  18:17:00  GMT.  The  storm  is  seen  to  be  small  in  ex¬ 
tent  (the  northwest  edge  is  not  observed  at  this  Cartesian  grid  level),  roughly 
2 

200  km  in  area.  It  is  composed  of  a  high  reflectivity  factor  cell  having  maximum 
value  near  40  dBZ.  The  storm  is  centered  approximately  104  km  west  and  70  km 
south  of  the  SPANDAR  radar.  Only  one  aircraft  penetration  was  recorded  this 
day.  It  covers  the  period  from  18:17:00  to  18:17:55  GMT,  and  the  time-adjusted 
track  is  shown  in  Figure  2.  The  observed  track  is  about  11  km  long  and  runs 
west  to  east.  It  begins  near  the  center  of  the  major  storm  cell  at  106  km  west 
and  68  km  south  of  the  radar,  about  3  km  north  of  the  highest  reflectivity  factor 
cell,  and  extends  to  the  radar  storm  boundary. 

The  sector  scan  radial  velocity  at  this  level  is  shown  in  Figure  3.  The  mag¬ 
nitude  range  is  roughly  -5  to  +2  m/s.  These  small  values  are  a  result  of  the  sec¬ 
tor  scan  being  oriented  nearly  orthogonally  to  the  environmental  wind.  It  is  seen 
that  the  wind  at  the  southwest  portion  of  the  storm  is  receding,  while  the  northeast 
portion  is  moving  towards  the  radar.  The  zero  velocity  line  running  southeast  to 
northwest  suggests  some  divergence  in  the  environmental  wind.  This  may 
indicate  that  the  storm  is  blocking  somewhat  the  onrushing  environmental  wind, 
causing  it  to  turn  northward  on  the  northeast  flank  of  the  storm.  The  major  cell 
penetrated  exhibits  a  component  towards  the  radar,  which  may  be  an  indication  of 
lower  level  air  rising  up  into  the  storm. 

The  radar  radial  velocity  time-history  plot  is  shown  in  Figure  4.  The  observed 
velocities,  with  the  small  positive  region  near  18:17:38,  agree  well  with  variation 
along  the  track  as  depicted  by  the  sector  scan  grid-point  velocities.  The  corre¬ 
sponding  tracking-gate  spectral  variance  estimates,  shown  in  Figure  5,  vary 

2 

from  about  2  to  6  (m/s)  .  The  variance  peaks  are  well  correlated  with  the  regions 
of  large  gradient  of  the  radial  velocity  (Figure  4).  The  somewhat  periodic  varia¬ 
tions  of  magnitude  in  the  velocity  and  variance  plots  exhibit  periods  in  the  range 
of  4  through  8  seconds.  With  an  aircraft  speed  near  250  mi/hr,  this  translates 
into  a  length  scale  of  450  to  900  m,  and  may  reflect  some  dominant  perturbation 

scale  regime  initiated  by  the  storm  structure. 
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The  turbulence  severity  estimates  (c  '  )  are  shown  in  Figur  e  The  5TF 
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derived  estimates  show  a  steady  inrrease  from  1.  5  (cm  '  )/s  at  18:17:25  GMT  to 

2  3 

a  fairly  constant  value  near  4  (cm  )  s  at  18:17:3  5  GMT  and  thereafter.  The 

2  3 

DSV  derived  estimates  start  near  5.5  (cm  )  s  and  naturally  vary  in  concert 
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with  the  spectrum  variance  (Figure  :>)  values  over  a  range  of  3.  5  to  4.  5  (cm  )  s. 
With  the  DSV  values  being  significantlv  greater-  than  the  ST  I-  estimates  prior  to 
18:17:35  GMT,  this  region  should  exhibit  rela'ivelv  low  shear,  have  most  turbu¬ 
lence  energy  mapped  into  the  DSV,  and  have  an  outer  scale  length  less  than  twice 
the  maximum  pulse  volume  dimension.  The  te-riod  stnt-t ine  at  18:17:35  GMT  shows 
quite  good  agreement  in  pattern  and  magnitude  lot  'he  t-..  q  estimates.  The  near 
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Figure  3.  Contours  of  Doppler  Radial  Velocity  in  m/s  for  1  July  1981  on 
Constant  Height  Surface  at  Aircraft  Storm  Penetration  Altitude  of 
4.57  km.  Contours  are  in  5-m/s  increments.  Time-adjusted  aircraft 
track  is  shown.  Distances  are  relative  to  SPANDAR  radar.  Positive 
values  represent  motion  away  from  the  radar.  Time  is  18:17:00  GT»'T 


1707  1715  1723  1731  1739  1747  1755 

TIME  IMIN/SECI 


Figure  4.  Time  History  Plot  of  Radar  Radial  Velocity  for  1  July  1981  as 
Measured  in  Tracking  Gate  During  Aircraft  Storm  Penetration.  Dashed 
lines  indicate  missing  data.  A  nine-point  block  filter  has  been  applied  to 
data 
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Figure  5.  Time  History  Plot  of  Doppler  Spectrum  Variance  (DSV)  for 
1  July  1981  as  Measured  in  Tracking  Gate  During  Storm  Penetration. 

A  nine-point  block  filter  has  been  applied  to  data.  Dash  line  repre¬ 
sents  missing  data 
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Figure  6.  Time  History  Plot  of  Turbulence  Severity  (e  >  )  as 
Estimated  by  Doppler  Spectrum  Variance  (DSV)  (Dash)  and 
Structure  Function  (Solid)  Techniques  for  1  July  1981 


equality  in  value  suggests  an  outer  scale  significantly  less  than  the  pulse  volume 
width  may  be  present.  However,  the  DSV  estimates  remain  roughly  constant 
(except  for  a  modulation),  while  the  STF  estimates  increase  to  equal  the  DSV 
values.  For  this  to  be  manifested  by  turbulence  alone  would  necessitate  e  to  vary 
as  an  unlikely  event.  A  more  reasonable  answer  is  that  the  fluctuation  of 

the  radial  velocity  (Figure  4)  is  more  representative  of  storm  wind  structure 
effects.  Then,  the  modulation  of  the  DSV  estimate  during  this  period  represents 
the  contribution  of  this  wind  shear  to  the  dominant  turbulence  component.  Thus, 
the  qualitative  conclusion  is  that  within  the  high  reflectivity  factor  storm  cell, 
little  radial  shear  is  present  and  most  turbulence  energy  is  mapped  into  the  DSV. 
As  the  storm  boundary  is  approached,  the  turbulence  severity  remains  relatively 
constant,  however,  storm-structure  wind  shear  increases.  The  effective  outer 
scale  length  is  less  than  1.8  km  throughout  and  the  DSV  appears  to  be  clearly 
superior  to  the  STF  method  used  here. 

5.2  3  July  1981 

On  3  July  1981  (day  184)  an  extensive  area  of  precipitation  moved  into  the 
region  from  the  southwest.  The  environmental  sounding,  Figure  7,  shows  both 
the  wind  direction  and  magnitude  to  be  relatively  constant  at  about  210°  and  12  to 
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15  m/s  from  4  km  to  9  km.  Figure  8  displays  the  storm  reflectivity  factor  at 

the  penetration  altitude  of  4.87  km  at  20:02:11  GMT.  The  region  of  interest  con- 
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sists  of  a  central  area  of  about  1800  km  and  is  located  near  80  km  west  and 
130  km  the  south  of  the  SPANDAR  radar.  The  area  contains  two  high  reflectivity 
factor  cells  with  maximum  values  near  45  dBZ.  The  time-adjusted  aircraft 
tracks  are  seen  to  extend  across  the  effective  storm  width  (20  to  20  dBZ  contour), 
and  pass  partially  through  the  two  storm  cells.  The  southwestern  cell  is  in  a  state 
of  decay.  At  the  time  of  20:15:32  GMT,  this  cell  has  essentially  dissolved  into 
the  background  with  a  maximum  reflectivity  factor  of  25  dBz. 
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I  iff.  ••  8.  (  'ini'iui  s  of  Storm  Reflectivity  Factor  at  Penetration  Altitude  of 

4.87  kn  for  i  Julv  1081  at  20:02:1  1  GMT.  Minimum  contour  value  is 


The  radial  velocity  values  in  Figure  9  range  from  about  0  to  -20  m/s.  Note 
that  the  radar  is  roughly  looking  along  the  wind  near  the  area  of  storm  penetration, 
with  the  aircraft  tracks  about  50  through  80°  off  the  wind.  The  pattern  of  radial 
velocities  show  increased  speed  along  the  southeast  and  western  boundary  of  the 
farther  cell,  with  a  minimum  value  on  its  upwind  (southwest)  side.  The  nearer 
cell  also  displays  increased  radial  wind  speed  along  its  eastern  boundary,  but 
decreased  speed  along  the  western  flank,  just  downwind  of  the  further  cell.  Al¬ 
though  the  wind  speed  maximum  along  the  eastern  edge  is  primarily  due  to  the 
radar  pointing  directly  along  the  wind  direction,  these  features  suggest  some 
blocking  effects  by  the  cells  are  occurring. 
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Figure  9.  Contours  of  Storm  Radial  Velocity  at  Penetration  Altitude  of  4.  87  km 
for  3  July  1981  at  20:02:11  GMT 


The  three  tracking-gate  radial  velocity  histories  are  shown  in  Figures  10 
through  12.  The  magnitude  ranges  from  about  -18  to  +1  m/s.  These  plots  show 
in  greater  detail  the  significant  modification  of  the  environmental  wind  on  the 
downwind  side  of  the  further  cell  (19:57:33  GMT),  and  the  receding  air  motion 
within  the  cells  themselves  (20:04:45  GMT).  These  data  agree  well  with  the 
larger  scale  grid-point  values  of  Figure  9. 

The  DSV  plots.  Figures  13  through  15,  show  the  spectrum  variance  to  be 
o 

generally  0.  5  to  2.  0  (m/s)  away  from  the  storm  cells,  but  to  have  rapid  and 
strong  variation  in  the  vicinity  of  the  cells.  For  example,  the  large  variation 
and  peak  magnitudes  from  19:57:53  to  19:58:25  GMT  of  run  1  coincide  with  pas¬ 
sage  through  the  near  storm  cell.  During  the  second  penetration  the  two  periods, 
20:04:29  to  20:04:55  GMT  and  20:05:30  to  20:05:50  GMT,  generally  coincide  with 
passage  through  the  two  cells.  The  periods  of  large  spectrum  variance  are  most 
often  not  associated  with  any  significant  gradients  of  radial  velocity.  This  sug¬ 
gests  that  pockets  of  strong  turbulence  rather  than  large-scale  wind  shear,  are 
most  often  the  cause  for  the  large  spectrum  variances  observed  here.  Last,  the 
time  histories  exhibit  fluctuations  with  per  iods  in  the  ranges  3  to  8,  13  to  16,  and 
24  through  35  sec,  or  r  oughly  700,  1500,  and  3500  m,  respectively.  The  short 
period  fluctuations  appear  as  a  r  ipple  on  the  large  period  fluctuations.  It  is  ob¬ 
served  that  the  energy  content  increases  dramatically  with  scale  size. 


Figure  10.  Time  History  Plot  of  Tracking  Gate  Radial  Velocity 
for  First  Penetration  at  4.87-km  .Altitude  for  3  July  1081 
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Figure  12.  Time  History  Plot  of  Tracking  Gate  Radial  Velocity 
for  Third  Recorded  Penetr  ation  at  4.  87-km  Altitude  on 
3  July  1981 


Figure  14.  Time  History  Plot  of  Tracking  Gate  Doppler 
Spectrum  Variance  (DSV)  for  Second  Recorded  Penetration 
at  4.  87-km  Altitude  on  3  July  1981 
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Figure  15.  Time  History  Plot  of  Tracking  Gate  Doppler 
Spectrum  Variance  (DSV)  for  Third  Recorded  Penetration 
at  4.  87-km  Altitude  on  3  July  1981 
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Figure  15.  Time  History  Plot  of  Tracking  Gate  Doppler 
Spectrum  Variance  (DSV)  for  Third  Recorded  Penetration 
at  4.  87-km  Altitude  on  3  July  1981  (Contd) 


The  turbulence  severity  estimates  are  displayed  in  Figures  16  through  18. 

The  maximum  pulse  volume  dimension  is  about  0.  9  km.  The  first  penetration 
track  is  about  55°  off  the  viewing  direction,  thus  the  apparent  shear  of  the 
radial  velocity  includes  contributions  from  both  longitudinal  and  transverse  com¬ 
ponents  of  storm  structure  wind  shear.  The  period  up  to  19:57:00  GMT  shows  the 
STF  and  DSV  estimates  to  be  small  and  nearly  equal  in  magnitude.  This  region, 
exhibiting  little  radial  shear  and  small  DSV,  suggests  that  turbulence  is  the  main 
contributor  to  spectrum  variance  and  the  outer  scale  length  is  much  smaller  than 
twice  the  maximum  pulse  volume  dimension,  thus  allowing  for  a  large  underesti¬ 
mate  in  the  DSV  value  to  approach  the  STF  value.  From  19:57:00  to  19:58:40  GMT 
the  DSV  estimates  are  considerably  larger  than  the  STF  values  and  the  two  patterns 
are  quite  similar,  suggesting  turbulence  i§  the  prime  contributor  to  spectrum 
variance  and  the  outer  scale  has  increased,  but  is  still  less  than  1.8  km.  The 
increasing  severity  from  19:57:00  to  19:58:00  GMT  suggests  relatively  constant 
with  increasing  turbulence  energy  (e)  as  the  aircraft  passes  from  downwind  of 
the  farther  storm  cell  (88  km,  122  km)  to  the  western  boundary  of  the  near  cell 
(78  km,  120  km).  The  DSV  turbulence  severity  within  the  near  storm  cell 
(19:58:00  to  19:59:00  GMT)  is  no  greater  than  that  found  at  the  cell  boundaries 
(25-dBZ  contour)  or  on  the  downwind  side  of  the  farther  cell  (19:57:30  GMT). 
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Figure  17.  Time  History  Plots  of  Turbulence  Severity  for  Second  Hun 
on  3  July  1981 
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Figure  18.  Time  History  Plots  of  Turbulence  Severity  for  Third  Run  on 
3  July  1981  (Contd) 


Passing  through  the  eastern  cell  boundary,  turbulence  severity  decreases  marked 

ly  with  the  outer  scale  length  being  less  than  1.  8  km.  The  period  near 

19:59:30  GMT  could,  once  again,  represent  a  situation  where  X  «  1.8  km  and 

o 

where  the  STF  estimate  is  strongly  biased  by  radial  shear. 

During  the  second  penetration  (Figure  17),  maximum  turbulence  severity  is 
found  at  the  eastern  edge  of  the  near  cell  and  center  of  the  farther  cell.  However, 
the  period  near  20:05:10  GMT,  corresponding  to  the  inter-cell  r  egion,  is  also 
highly  turbulent.  The  patterns  of  the  two  estimates  are  quite  dissimilar  through¬ 
out  the  near  cell,  while  in-phase  through  the  more  distant  cell.  Beyond  the 
farther  cell,  the  two  estimates  are  roughly  equal  in  magnitude  but  differ  in 
pattern.  These  observations  suggest  that  regio  is  of  str  ong  storm  wind  shear 
within  the  near  cell  are  biasing  the  turbulence  measurements,  in  particular,  the 
STF  estimates  (for  example,  20:04:20  GMT).  Little  storm  structure  shear  and 
strong  turbulence  are  expected  in  the  farther-  cell.  The  outer  scale  lengths  are 
expected  to  be  less  than  1.  8  km  thr  oughout.  I  pwind  of  the  farther  cell  shear  of 
the  storm  wind  field  may  be  contributing  significant  ly  to  the  turbulence  estimate's. 

The  third  penetration  (Figure  18),  in  reasonably  close  proximity  to  the  first 
penetration,  exhibits  behavior  consistent  with  earlier  observations.  These  are 
moderate  turbulence  downwind  of  the  farther  cell  (20:1  1:00  C.M  T),  strong  turbu¬ 
lence  at  the  western  edge  (20:11:40  GMT),  and  within  (20:11:57  GMT)  the  near  cell 


The  strongly  turbulent  period  at  12:21:00  GMT  corresponds  to  the  nondescript 
25  dBZ  region  east  of  the  near  cell. 

The  period  near  20:12:09  GMT  shows  a  strong  drop  in  the  DSV  estimate 
coincident  with  a  strong  rise  in  the  STF  value.  Figure  12  shows  a  region  of 
moderate  shear  of  the  radial  velocity  at  this  time.  This  may  be  a  period  where 
turbulence  strength  (e)  decreases  coincidently  with  a  strong  increase  in  storm 
structure  wind  shear.  Noting  that  the  STF  estimates  during  this  time  generally 
include  part  or  all  of  the  radial  velocity  data  between  20:12:00  to  20:12:21  GMT, 
this  may  be  a  case  where  storm  structure  is  strongly  influencing  the  STF  estimates 
and  may  be  presenting  an  inaccurate  picture  of  the  turbulence  field  at  this  loca¬ 
tion.  This  is  particularly  so  since  the  DSV  estimate  does  not  respond  in  kind. 

These  data  show  a  complex  distribution  of  moderate  to  heavy  turbulence. 

The  turbulence  severity  is  high,  not  only  within  the  storm  cells,  but  also  at  the 
storm  cell  boundaries,  in  between  the  two  cells,  and  in  nondescript  2 5 -dBZ 
regions  outside  and  downstream  of  the  storm  cells.  The  DSV  estimates  are 
generally  much  greater  than  the  STF  values,  indicating  A  less  than  twice  the 
maximum  pulse  volume  dimension  of  0.9  km.  The  DSV  estimates  provide  a 
clearer  picture  of  the  turbulence  strength,  while  the  STF  estimates  appear-  to  be 
unreliable  due  to  inclusion  of  storm  structure  wind  shear. 

5.3  16  July  1981 

On  16  July  1981  (day  197)  a  broad  area  of  precipitation  moved  into  the  ob¬ 
servational  region  from  the  west.  The  sounding  data  is  shown  in  Figure  19.  The 
environmental  wind  speed  varies  from  near  0  m/s  at  the  sur  face  to  10  m's  near 
3.  5  km,  and  remains  at  about  12  m/s  up  to  6  km.  At  the  surface  the  wind  direc¬ 
tion  is  roughly  200°,  changes  to  270°  at  2-km  height,  and  remains  nearly  constant 
at  280  through  310  above.  Figure  20  shows  a  reflectivity  factor  plot  at  the  F106B 
penetration  height  of  4.27  km  for  the  time  18:32:37  GMT.  The  two  adjusted  air¬ 
craft  tracks  are  also  shown.  The  area  of  interest  is  a  high  reflectivity  factor- 
cell  having  a  maximum  value  near  40  dBZ  on  the  western  edge  of  the  storm  com¬ 
plex.  It  is  located  at  113  km  west  and  130  km  south  of  the  SRANDAK  radar  and 
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comprises  an  area  of  about  400  km  .  The  first  recorded  aircraft  track  begins 
near  the  center  of  the  storm  cell  and  extends  southward  to  the  storm  boundary. 

The  second  track  flirts  with  the  western  edge  of  the  storm. 

The  storm  radial  velocity  structure  is  shown  in  Figure  21.  The  velocity 
magnitudes  are  relatively  stable  near  5  m  s,  except  for  a  small  region  where  the 
oncoming  radial  wind  is  greater  than  10  in  s,  on  the  southern  edge  of  the  cell. 

With  the  environmental  wind  from  the  northwest  the  two  aircraft  tracks  and  the 
radar  viewing  direction  are  essentially  orthogonal  to  it,  thus  the  small  radial 
velocities  observed.  A  later  scan  at  1 1 ' : 0 0 : ">  12  GMT,  showing  similar  si/e,  shape. 
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Figure  20.  Contours  of  Storm  Reflectivity  Factor  at  Penetration 
Altitude  of  3.  5  km  for  16  July  1981  at  18:32:37  GMT.  Minimum 
contour  value  is  15  dBZ 
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Figure  24.  Time  History  Plot  of  Tracking  Gate  Doppler  Spec¬ 
trum  Variance  (DSV)  for  First  Recorded  Penetration  at  3.  5 -km 
Altitude  on  16  July  1981 

Figure  26  displays  the  turbulence  severity  for  the  second  storm  boundary 
penetration.  The  block-averaged  variance  data  obtained  during  the  first  pene¬ 
tration  contained  spurious  values  and  were  not  used  in  this  analysis.  The  DSV 
derived  estimates  exhibit  a  simple,  slowly  decreasing,  turbulence  severity.  The 
STF  derived  values,  on  the  other  hand,  show  a  contradictory  increasing  magnitude 
over  the  first  half  of  the  history,  while  good  agreement  exists  over  the  latter  half. 
The  ever-increasing  STF  values  are  believed  to  result  in  part  from  the  STF 
analysis  itself.  Nonetheless,  the  features  suggest  mild  turbulence  with  increasing 
wind  shear  during  the  first  half  of  the  penetrations,  with  mild  turbulence  and  little 
shear  afterwards.  The  data  suggest  the  outer  scale  length  is  probably  less  than 
twice  the  maximum  pulse  volume  dimension  of  1.3  km. 

5.4  17  July  1981 

The  final  day  of  observations  was  17  July  1981  (day  198).  The  sounding  data 
are  shown  in  Figure  27.  The  environmental  wind  speed  increases  slowly,  from 
about  2  m/s  near  1  km  height  to  about  15  m  s  at  5  km  height.  It  then  is  nearly 
constant  at  14  m/s  up  to  7  km,  where  it  then  gradually  increases  to  28  m/s  at 
11  km.  The  wind  direction  varies  from  nearly  130°  at  the  surface  to  300°  at 
2.0  km,  and  then  remains  constant  above.  The  storm  reflectivity  factor  is 
displayed  in  Figure  28  for  the  aircraft  penetration  altitude  of  3.35  km.  It  is  a 
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Figure  27.  Sounding  Data  for  17  July  1981  at  1200  GRIT 


small  complex,  with  the  area  of  interest  located  about  80  km  due  east  of  the 
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SPANDAR  radar.  This  nearest  storm  is  only  about  100  km  and  is  composed  of 
a  single  high  reflectivity  factor  cell  having  magnitude  near  50  dBZ.  The  plot  time 
is  18:35:00  GRIT,  The  three  adjusted  aircraft  tracks  show  the  aircraft  penetrated 
the  southern  and  northern  boundaries  of  the  storm  in  regions  where  the  reflectivity 
factor  is  25  dBZ  or  less. 

The  corresponding  radial  velocity  plot  is  shown  in  Figure  29.  This  velocity 
structure  is  undeniably  the  most  complex  of  the  four  days  presented.  The  general 
trend  of  increasing  magnitude  towards  the  southern  boundary  reflects  the  radar 
scanning  more  closely  along  the  environmental  wind  direction.  The  region  of 
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Figure  28.  Contours  of  Storm  Reflectivity 
Factor  at  Penetration  Altitude  of  3.35  km 
for  17  July  1981  at  18:35:00  GMT.  Mini¬ 
mum  contour  value  is  15  dBZ 


slight  increase  in  radial  wind  speed,  having  the  appearance  of  the  horseshoe¬ 
shaped  5  to  10  m/s  region  upwind,  adjacent  to,  and  downwind  of  the  cell,  do 
however,  suggest  some  divergence  of  the  environmental  wind  as  it  encounters  the 
storm.  The  penetration  tracks  are  roughly  20°  off  the  wind  direction. 

The  time  histories  of  tracking-gate  Doppler-  mean  velocity  and  variance  are 
shown  in  Figures  30  through  32  and  Figures  33  through  35,  respectively.  The 
radial  velocity  histories  are  in  general  agreement  with  the  grid-point  values. 
However,  the  first  track  shows  a  reversal  in  wind  dir  ection  near  9  km  south  and 
87  km  west,  between  the  major  storm  and  the  smaller  storm  further  west.  It  is 
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Figure  29.  Contours  of  Storm  Radial  Velocity  at 
Penetration  Altitude  of  3.35  km  for  17  July  1981 
at  18:35:00  GMT. 


believed  that  the  pulse-pair  velocities  here  are  biased  due  to  low  returned  power 
and  are  considered  suspect.  The  same  is  true  for  the  region  of  negative  values 
located  7  km  north  of  this  location.  In  general,  the  data  show  a  range  of  velocities 
from  -2  to  9  m/s.  The  periods  of  missing  data  in  run  3  indicate  the  track  was 
skirting  the  radar  detectable  northern  boundary.  Tracks  that  appear  outside  the 
plotted  storm  boundaries  were  generally  in  regions  having  reflectivity  factor  of 
5  to  15  dBZ. 
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Figure  31.  Time  History  Plot  of  Tracking  Gate  Radial  Velocity 
for  Second  Recorded  Penetration  at  3.35-km  Altitude  for 
17  July  1981 
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Figure  32.  Time  History  Plot  of  Tracking  Gate  Radial  Velocity 
for  Third  Recorded  Penetration  at  3,35-km  Altitude  for 
17  July  1981 
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Figure  32.  Time  History  Plot  of  Tracking  Gate  Radial  Velocity 
for  Third  Recorded  Penetration  at  3.35-km  Altitude  for 
17  July  1981  (Contd) 
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Figure  33.  Time  History  Plot  of  Tracking  Gate  Doppler  Spec¬ 
trum  Variance  (DSV)  for  First  Recorded  Penetration  at  3.35-km 
Altitude  on  17  July  1981 
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Figure  33.  Time  History  Plot  of  Tracking  Gate  Doppler  Spec- 
,  trum  Variance  (DSV)  for  First  Recorded  Penetration  at  3.35-km 

*  Altitude  on  17  July  1981  (Contd) 


Figure  34.  Time  History  Plot  of  Tracking  Gate  Doppler  Spec¬ 
trum  Variance  (DSV)  for  Second  Recorded  Penetration  at 
3.35-km  Altitude  on  17  July  1981 
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Figure  34.  Time  History  Plot  of  Tracking  Gate  Doppler  Spec 
trum  Variance  (DSV)  for  Second  Recorded  Penetration  at 
3.35-km  Altitude  on  17  July  1981  (Contd) 


Figure  35.  Time  History  Plot  of  Tracking  Gate  Doppler  Spec 
trum  Variance  for  Third  Recorded  Penetration  at  3.35-km 
Altitude  on  17  July  1981 
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The  variance  time  histories  display  significant  fluctuations,  particularly  dur¬ 
ing  the  latter  two  penetrations.  This  is  in  agreement  with  the  greater  fluctuation 
observed  in  velocity.  One  must  keep  in  mind,  however,  that  because  of  the  low 
signal-to-noise  ratio  occurring  during  these  two  runs,  larger  scatter  about  the 
true  Doppler  spectrum  mean  and  variance  are  to  be  expected.  As  in  previous 
data,  regions  of  large  spectrum  variance  are  only  occasionally  correlated  with 
strong  gradients  of  the  radial  velocity.  Remaining  aware  of  the  potential  problems 
from  low  signal  power,  these  data  suggest,  nonetheless,  that  the  large  variance 
values  observed  result  more  from  strong  turbulence  than  from  shear  of  the  large- 
scale  radial  wind.  Last,  observations  of  the  histories  suggest  that  the  large  var¬ 
iance  values  observed  exhibit  fluctuations  with  preferred  periods  in  the  range  of 
9  to  12  sec,  which  translate  into  length  scales  of  about  1. 0  to  1.5  km. 

Figures  36  through  38  display  the  turbulence  severity  estimates  for  the  three 
recorded  penetrations  for  this  day.  The  pulse  volume  width  is  about  0.  6  km.  The 
first  pass  shows  reasonably  good  agreement  between  the  two  estimates  during  the 
first  half  of  the  track.  Through  the  second  half  the  DSV  estimates  increase 
dramatically,  while  the  STF  estimates  remain  nearly  constant.  The  STF  data 
suggests  that  wind  shear  was  not  a  significant  contributor  to  DSV.  However, 
since  the  aircraft  flew  essentially  along  the  radar  viewing  direction  for  these 
flights,  no  true  measure  of  the  transverse  storm  wind  shear  is  available.  These 
data  may  also  indicate  that  the  significant  broadening  of  the  Doppler  spectra  is 
due  to  low  signal  power,  however ,  one  would  also  expect  the  STF  estimate  to 
increase  if  this  were  the  case.  Thus,  the  data  may  suggest  a  small  effective 
outer  scale  length,  generally  increasing  in  size,  but  still  less  than  1.2  km  as  the 
plane  moves  west  to  east. 

The  second  penetration  exhibits  similar  behavior  when  the  DSV  estimate 
remains  relatively  constant  while  the  STF  values  reach  a  maximum  near 
18:39:49  GMT.  However,  the  strong  radial  shear  in  the  radial  wind,  as  is  shown 
by  Figure  31,  suggests  that  the  STF  estimate  is  biased  by  storm-structure  wind 
shear. 
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Figure  3(1.  Time  History  Plots  of  Turbulence  Severity  for  First  Hun 

on  17  July  1981 


Figure  37.  Time  History  Plots  of  Turbulence  Severity  for  Second  Hun 

on  17  July  1981 


6.  CONCLUSIONS 


The  data  presented  here  describe  the  basic  observations  made  during  the 
1981  Joint  Agency  Turbulence  field  experiment.  A  total  of  four  aircraft  storm- 
penetration  periods  were  recorded,  and  resulted  in  nine  spatially  and  temporally 
coincident  tracking  Doppler  radar  and  in-situ  aircraft  data  sets.  A  variety  of 
environments  ranging  from  nearly  uniform,  mildly  turbulent  wind  fields,  to 
highly  turbulent  and  complex  fields  were  observed.  In  most  cases  the  DSV  fluc¬ 
tuations  observed  in  the  tracking  gate  data  are  generally  not  well  correlated  with 
strong  gradients  of  the  radial  wind  field.  This  suggests  that  turbulence,  rather 
than  wind  shear  is  primarily  responsible  for  the  large  spectrum  variance  observed. 
The  time  history  data  also  show  persistent  fluctuations  with  a  range  of  scales  from 
200  m  to  3.  5  km,  with  the  energy  content  increasing  dramatically  with  increasing 
scale  length.  Estimates  of  turbulence  severity  show  significant  turbulence  dis¬ 
tributed  throughout  the  storm  regions.  Heavy  turbulence  is  often  found  not  only 
within  high  reflectivity  factor  storm  cells,  but  also  in  between  cells,  near  storm 
boundaries,  and  in  innocuous -appearing  low  reflectivity  factor  regions.  Observa¬ 
tions  suggest  that  the  turbulence  fields  are  generally  characterized  by  scales  of 
motion  less  than  about  2  km  in  length.  At  long  range,  DSV  methods  are  definitely 
superior  to  structure  function  techniques  that  do  not  account  for  beam  filtering 
effects.  Storm -structure  wind  shear  may  strongly  bias  the  structure  function 
estimates,  but  generally  do  not  significantly  influence  the  variance  methods. 

Similar  analyses  will  be  performed  on  the  aircraft  data,  and  these  data  will  serve 
as  'ground  truth'  for  the  radar  data.  These  efforts  will  be  used  to  determine  the 
energy  and  spatial  characteristics  of  the  turbulent  fields,  and  allow  for  adjust¬ 
ments  in  the  radar  methodology.  These  turbulence  analyses  will  appear  in  a 
later  report. 
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